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Abstract: We point out that in some situations it is possible to use matrix model tech- 
niques a la Dijkgraaf-Vafa to perturbatively compute D-brane instanton effects. This pro- 
vides an explanation in terms of stringy instantons of the results in hep-th/0311181. We 
check this proposal in some simple scenarios. We point out some interesting consequences 
of this observation, such as the fact that it gives a perturbative way of computing stringy 
multi- instanton effects. It also provides a further interpretation of D-brane instantons as 
residual instantons of higgsed supergroups. 
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1. Introduction 

In a series of remarkable papers [jy, [2|, g], Dijkgraaf and Vafa have argued that the exact 
low energy superpotential for a class of gauge theories can be computed by taking the 
planar limit of a matrix model. This matrix model is determined in a simple way by the 
deconstruction 0] of the original supersymmetric theory. It has a set of (bosonic) matrix 
valued fields in one to one correspondence with the superfields of the original theory, and a 
potential given by the superpotential of the original theory. For completeness, we include 
a short review of the Dijkgraaf-Vafa results in Section |2.l| . 

Remarkably, even if the effective superpotential is obtained from the planar limit of 
the matrix model, it is the low energy superpotential for finite iV in the original theory 
(with N not necessarily large). A particularly interesting case is that in which N is small. 
This case led originally to some puzzles, since for some particular cases the results obtained 
from the matrix model and the pure gauge theory seemed to disagree. The situation was 
clarified in [||, |6[, where it was argued that the matrix model results might indeed differ 
from the field theory expectations for gauge factors of low rank, but the disagreement could 
be attributed to the choice of a UV definition of the gauge theory: the matrix model gives 
results for the definition of the gauge theory corresponding to the particular geometric 
embedding of the gauge theory in string theory chosen by Dijkgraaf and Vafa. 

In H |6| the UV completion of the gauge theory was defined by its embedding in a 
(higgsed) supergroup of high enough rank. At the level of F-terms this definition of the 
gauge theory coincides with the usual gauge theory one, but it has the advantage that it 
is well defined in the UV. It is also very natural from the stringy point of view, where it 
can be understood as adding to the original DV embedding of the gauge theory an infinite 
set of brane-antibrane pairs with trivial K-theory charge. 
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The purpose of this note is to point out, in the simplest possible set of examples, that 
it is also possible to understand this discrepancy between gauge theory and matrix models 
as coming from stringy D-brane instanton effects |?], ||, [|. Our claim is that the matrix 
model computation includes the contribution of D-brane instantons to the low energy 
superpotential of the string theory realization of the gauge theory. 

More precisely, we claim that it is possible to obtain the contributions of D-brane 
instantons to the effective superpotential by studying the planar limit of the associated 
matrix model. In the conclusions we include some preliminary remarks about the interpre- 
tation of the non-planar contributions, and we also make some comments on how could it 
be possible to generalize our results to some geometries beyond the Dijkgraaf-Vafa setup. 

An interesting consequence of this result is that it gives a relatively easy and systematic 
way of computing perturbatively stringy multi-instanton effects, at least when a matrix 
model description of the system is available. These multi-instanton effects are typically 
rather involved to study using stringy instanton calculus, so the perturbative approach 
based on the matrix model will be a welcome tool. 

We can find some other interesting implications of the result. There are by now differ- 
ent ways of computing exotic D-brane instanton effects by modifying the system such that 
it admits a gauge theory interpretation while still allowing to extract the relevant informa- 
tion. Let us mention as examples [1C, 11, 12], in which Seiberg duality (as in the duality 



cascade) is used to turn D-brane instantons into conventional gauge non-perturbative ef- 



fects. Another example is [13], where Seiberg duality is combined with higgsing. 



From this perspective, the matrix model viewpoint offers yet another embedding of D- 
brane instanton effects into gauge theory. It is given in terms of the supergroup construction 
of []|], which essentially consists of adding infinite brane-antibrane pairs to the system. In 
this setup stringy D-brane instanton effects are identified with residual instanton effects of 
the higgsed supergroup. 

Some of the constructions above also admit a nice interpretation in terms of matrix 
models. For example, as explained in [||, Seiberg duality in the case of ADE quiver theories 
can be understood as Weyl reflections of the quiver diagram, and the matrix model results 
are essentially insensitive to this Weyl reflection. Similarly, as described in Q, the F-terms 
computed by the matrix model are insensitive to where we are in the duality cascade, so 
we can do our calculations either "up the cascade", where everything can be formulated 
in terms of gauge theory, or at the bottom of the cascade, where we have exotic D-brane 
instanton effects with no gauge theory interpretation. Let us mention that [Q also provides 
a nice intermediate viewpoint between the brane-antibrane configuration and the matrix 
model in the form of the stiperaiatrix model description, which turns out to be rather useful 
for understanding the physics of these duality cascades in the context of matrix models. 

We would also like to point out that very similar considerations to some of the ones we 



discuss here have already appeared in [14]. In that paper, the geometric transition is used 



to compute non-perturbative effects due to stringy D-brane instantons. The matrix model 
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and the geometric transition approaches are very closely related, and they can be shown 
to give the same result in the planar limit, which determines the superpotential. This 
is, nevertheless, a large iV effect 1 , and we expect that we have to use the matrix model 
for computing higher derivative interactions in the effective superpotential. If we restrict 
ourselves to the superpotential, and thus the planar limit, the matrix model provides a 
different but equivalent viewpoint of the discussion in [14]. 



We feel that the previous reasons make it worth pointing out the explicit connection 
between D-brane instantons and matrix models. It is, in any case, a subtle and interesting 
connection between two rich fields of study. 

Let us mention that in this note we will be checking the proposed connection only in the 
simplest possible scenario: we will be computing one-instanton superpotential contributions 
for quiver gauge theories admitting a matrix model description a la Dijkgraaf and Vafa. 
This leaves many issues to be studied, and we hope to come back to them in future work. 

The remainder of this note is organized as follows. We start in Section ^ by summariz- 
ing the main points of the Dijkgraaf- Vafa correspondence and D-brane instanton calculus. 
In Section || we see how the D-brane instanton and matrix model results agree for the 
simplest one-instanton sector in conifold geometries. In Section || we study a more inter- 
esting and involved example, the abelian orbifold of the conifold. Section |5] contains our 
conclusions and some discussion of interesting open problems. 

2. Background material 

2.1 The Dijkgraaf- Vafa correspondence 

Let us start by reviewing the main results of Dijkgraaf and Vafa jl], ^, ||. Since we will be 
dealing with stringy instantons we concentrate on the stringy derivation and interpretation 
of the results, although the same results can also be derived in pure field theory Jl|, |l|, . 



Gauge theory 

The configurations we will be concerned with are stringy embeddings of N = 1 gauge 
theories with a mass gap, the lightest states being glueball superfields 2 . The canonical 
example is a N = 1 theory obtained from J\f = 2 U(N) SYM with no flavors and a 
potential for the adjoint superfield: 

d 2 9TrW($) (2.1) 

This superpotential breaks down the supersymmetry to Af = 1. We take W(z) to be a 
polynomial function of degree n + 1. We will denote the dynamical scale of this theory by 
A. 



1 We hope that the terminology we use does not lead to confusion. The N here refers to the rank of the 
matrices in the matrix model, and has nothing to do with the ranks of the physical gauge group. We review 
this point in Section 



2.1 



2 The existence of the mass gap, and glueballs being the correct low energy degrees of freedom, are 
statements which are very hard to show rigorously. Here we just assume them to hold. 
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The classical vacua of these theory are obtained by extremizing the superpotential W. 
This superpotential has n extrema ai, let us indicate this by writing the derivative of W 
as W'(z) = g]Xi=i(. z ~ a i)- A classical solution of Tr W'{<&) = is obtained by considering 
a partition of the iV eigenvalues of $ into the n extrema: 

Ai \ 



\ 



Around any of these vacua $ is massive so it can be integrated out of the low energy theory, 
which is then given by Af = 1 SYM with the Higgsed gauge group YYl=i U(Ni). The Ni 
here count how many eigenvalues of $ are equal to a^. It is conventional and convenient 
to denote this Higgsing process U(N) — > Yl™—i U(Ni). 

This theory confines, and it is believed that the low energy degrees of freedom of this 
theory after confinement, similarly to the case of pure SYM, are given by the glueball 
superfields Si, one for each of the gauge factors: 

1 



(2.3) 



with W a the spinorial chiral superfield appearing in the gauge kinetic action, having lowest 
component the gluino X a . As explained by Veneziano and Yankielowicz [18|, for pure SYM 
the low energy dynamics can be obtained from an effective superpotential involving the 
glueball superfield: 

W eff (S) = Slog (A 3N /S N ) + NS (2.4) 

where A denotes the dynamical scale of the SU(N) theory. 

In the case of the theory we are dealing with, the effective superpotential will be of 
the form |nj|: 



Weff(S)=J2 



Si log 



A 



3Ni 



s 



N, 



+ + Ni 



QMS) 

dS,; 



(2.5) 



where J-q is a perturbative series on the glueballs S%, and Aj is the dynamical scale of the 
factor SU(Ni) after integrating out the massive matter. The usual matching relations lets 
us express Aj in terms of the dynamical scale A of the original theory, the derivative of the 
superpotential W'(z), and the residual ranks Nf. 



A 



3JVi 



A 2N gNi 



ai) 



Ni-2Ni 



(2.6) 



Similar matching relations for the related SO(N) — ► SO(Nq) x Y\ i U(Ni) and Sp(N) — > 
Sp(No) x Y\ i U(Ni) cases can be found in @. 

What Dijkgraaf and Vafa showed is that J-q can be computed perturbatively from a 
related matrix model. We describe the main points of their argument below, but let us 
first describe how we will embed these gauge theories into string theory. 
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Geometry 

A type IIB configuration realizing the gauge theory we are studying was described in [19] 



It is given by Minkowski spacetime times a local Calabi-Yau defined by the resolution of: 

x 2 + y 2 + w 2 + {W'{z)f = (2.7) 



with W{z) the same function as in (2.1), and W'(z) its first derivative with respect to z. 



The extrema of W(z) are located at z = ai, so we can rewrite ( |2.7|) more clearly as: 

n 

+ y 2 + w 2 + g 2 H(z - en) 2 = 0. (2.8) 



x 2 



i=l 

It is easy to see from here that when z is close to aj (we take the a% to be well separated 
for simplicity), the equation (|2.Sf ) reduces to the one for the conifold. We take the conifold 
to be resolved 3 , so the geometry is no longer singular, but instead has an isolated 2- 
cycle at each of the degenerations a». The size of the resolved 2-cycles is determined by 
the complexified gauge coupling constant for the theory we are engineering, more on this 
below. 

Let us consider this background in the presence of N D5 branes, which we take to wrap 
Minkowski and two internal directions. The low energy dynamics of such a system is U(N) 
with the superpotential ( |2,1| ) for the adjoints. Classical vacua of the system correspond to 
distributing the N D-branes in different ways in the n different 2-spheres of the resolved 
conifolds. When the N branes get separated into the different cycles, some open strings 
get massive, and the gauge group on the branes factorizes. Furthermore the adjoint chiral 
multiplet living on the brane, and which parametrizes motions of the brane, gets a positive 
mass since the resolution two-sphere is a local volume minimum in its homology class. 

Let us briefly mention that all the resolved two spheres belong to the same homology 
class. One simple way to see this in physics terms is that we can tunnel from U (iVi) x U (N2) 
to U(Ni — 1) x U(N2 + 1) in the gauge theory just by smoothly taking one eigenvalue of 
<I> from ai to a%. For such a tunneling to be possible in string theory the total homology 
class should stay the same, the a\ — > a% interpolation describes the 3-chain that makes the 
two cycles homologous. Such a tunneling is described by instantons of the matrix model 
|l|, J2l| , which are related to domain walls in the effective theory interpolating between the 
two vacua, and coming from branes wrapping the 3-chain. 

The bare holomorphic coupling constant r of the Yang-Mills theory living on the brane 
is given by the complexified volume of the two cycle S wrapped by the brane: 

4vrz 1 f , „. , 

{J + iB). (2.9) 



2vr g\ M 2mg s J s 

We will also be interested in engineering Sp(N) theories with a superpotential for the 
adjoint. They can be obtained by adding a suitable orientifold to the geometry above. The 
breaking pattern in this case is Sp(N) — > Sp(Ni) x U(N2) x . . . x U(N n ), but otherwise 
everything is very similar to the U(N) case ||. 

3 Resolution means, roughly, that we substitute the singularity of the conifold by a finite size S 2 in such 
a way that the resulting space is smooth. See for example po| for a more detailed description. 
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The low energy degrees of freedom of the string theory 

We see that up to now the tree level dynamics of the stringy system and the gauge theory 
are very similar, and can be matched quite easily. Nevertheless, we have not yet described 
what plays the role of the all-important glueball superfield in the string side. This is a 
subtle problem, and central to our discussion. 

It turns out that in order to account for all the low energy dynamics of string theory 
one needs to include more "glueball" superfields in the string theory side than those that 
would be necessary just from the gauge theory point of view. In fact, one way of reading 
the results of this note is that these extra degrees of freedom in the string side are those 
necessary in order to account for the stringy instanton effects, which are not present in 
gauge theory. 

These extra glueballs must be included for the cases in which there is a single brane 
wrapping a S 2 , giving rise to a U(l) gauge group, and in the case on which the S 2 is 
wrapped by no branes, but there is an orientifold such that the gauge group on the cycle, 
had there been N branes wrapped on it, would have been Sp(N). We will refer to this last 
case as the Sp(0) case. In gauge theory, none of this two gauge groups would give rise to a 
glueball and a Veneziano-Yankielowicz superpotential describing its dynamics. These are 
also the only cases on which one has to include these extra degrees of freedom in order to 
match with the superpotential due to one stringy instanton effects. Other candidates one 
could think of such as 50(0), 50(2) or U(0) do not require extra glueballs for describing 
the stringy dynamics. The justification of this result from our point of view is that these 
are the D-brane instantons that contribute to the superpotential (see Section |2^ below), 
so we need to make the corresponding glueballs dynamic in order to reproduce the known 
stringy instanton effects. 

Note that the prescription given in for when to include the extra stringy glueballs, 
which was obtained via the geometric transition, coincides with our prescription that the 
matrix model should compute D-brane instanton effects, which is a successful check of our 
proposal. 

We should point out that the previous discussion is the correct one when we are 
restricting ourselves to superpotential contributions in the single instanton case (as we are 
doing in this note), but we expect that it might be necessary to consider dynamical at 
least the U(0) glueballs in the description of more complicated processes. The reason is 
the following. Consider instantons in a U(0) gauge group (these are commonly called U{\) 
instantons in the D-brane instanton literature). They have four neutral zero modes, and 
thus they do not contribute to the superpotential by themselves, in agreement with the 
fact that we are considering the corresponding glueball non-dynamical. Nevertheless, as 
discussed in multi-instanton processes involving a Sp(0) and a U(0) instanton can 
contribute to the superpotential. We cannot reproduce these effects in the matrix model 
side if we just set the U (0) glueball to zero from the beginning. One way out is to include 
the U (0) glueball in the game, and just set to zero the Veneziano-Yankielowicz part of the 
superpotential. Namely, we take the following formula to hold for any N > 0: 

W U(N) = 5 log (Al N /S N ) +NS- 2-kItS (2.10) 
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Here Ao is the cutoff for the theory, which combines with the bare holomorphic coupling 
r = J 4 £^- in order to give the usual holomorphic scale A of the gauge theory: 

A 3 = A 3 , e ~ 2niT/N (2.11) 

reproducing in this way the Veneziano-Yankielowicz superpotential fl2.4Q . 

Note that the first two terms in ( |2.10| ) vanish when N = 0. For the one-instanton 
superpotential, such a prescription dynamically sets 5 = due to the effect of an extra 
light hypermultiplet appearing close to S = ||, but for multi-instanton processes there 
could be nontrivial contributions. 

This prescription for the U(0) glueball has a chance of matching the multi-instanton 
results, and also fits well with the heuristic idea that dynamical glueballs are associated 
with instanton effects in the corresponding node. We hope to come back to this interesting 
question in the future. In any case, we will not go beyond the one-instanton sector in this 
note, so this subtlety will not be important and we refrain from further discussion of this 
point. 

Matrix models as the theory on the branes 

Let us now come back to the connection between the string theory and the matrix model, 
following |2|, H|. Since we are interested in computing F-terms, we can go to the topo- 
logical string. As described above, we have N D5 branes in a geometry defined by the 
function W(z). The topological theory living on the worldvolume of the D5 branes is 
the reduction of three complex dimensional holomorphic Chern-Simons theory |2^] to the 
two dimensional worldvolume of the brane. Partially solving the equations of motion for 
this theory tells us that the F-term physics is captured by the partition function of the 
holomorphic matrix model with action: 

2= ———— 1 TT/nr .. exp ( -— / d$TrW($)) (2.12) 

Vol(C7(Mi) x • • • x U(M n )) y \g s J V ') y 1 

This formula requires some explanation. Here is a M x M matrix (M is not related to 
the number of branes N, its relevance to the physical theory will be explained shortly), and 
g s is a overall scale unrelated to the coupling constant of the physical theory. As above, the 
classical vacua of this theory will be determined by extremizing VK($), and distributing 
the M vacua into the n extrema. We take M = M\ + M2 + . . . + M n . The prefactor takes 
into account the volume of the unbroken gauge group in the chosen vacuum. 

For the particular case of the superpotential, which will be the main focus of this note, 
it can be further argued [^] that only the planar diagrams contribute 4 . This means that 
the required information can be obtained in the large M limit of the matrix model. As 
usual when taking large M limits, we will keep the 't Hooft coupling Si = g s Mi fixed. This 
has two desirable consequences. First it tells us that we can solve the matrix model in 

4 The argument boils down to the fact that the planar diagrams are the ones without momentum inser- 
tions. It also implies that the higher F-terms, which do involve momentum insertions, are computed by the 
non-planar part of the matrix model. We will come back to this observation in the conclusions. 
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the saddle point approximation, which is easier; and second, one can apply the arguments 



of 24, [25|] to argue that the result can be computed via a geometric transition to a flux 



background. 

Let us now describe how the prescription goes for connecting the matrix model ampli- 
tudes to F-terms in the physical string. The 't Hooft parameters Si of the matrix model 
get identified with the glueballs Si of the physical string. In this way, the perturbative 
expansion of the matrix model is an expansion in powers of the glueballs. The conjecture 
is that the low energy superpotential of the string theory is given by: 



t=l 



oSi 



(2.13) 



where Ni denotes the rank of the gauge factor U(Ni) and t% is the bare holomorphic 
gauge coupling on the branes wrapping the cycle which is determined by the volume of 
the wrapped two cycle as in eq. ( |2.9D 5 . Zq denotes the planar part of the matrix model 
partition function ( 2.12[) . 



Expressions (|2.13[) and ( |2.5| ) match nicely once we take into account the following. 



As explained in [26] from the mirror Chern-Simons perspective, the answer for the free 
planar part of the matrix model, once we take the volume of the residual gauge group into 
account, gives the Veneziano Yankielowicz superpotential. We can then expand the planar 
partition function into a free part and a correction: 

exp ( -— f d$ Tr ^)\ + T (2.14) 



Vol(Z7(M 1 ) x • • • x U{M n )) \ g s J 2 

The first term gives rise to the Veneziano- Yankielowicz term in the superpotential (notice 
also that it is the partition function for the conifold), and J-q is a polynomial in Si'. 

■^0 = E 'm...^S'---^ (2-15) 

11. ..i n 

This series expansion in S{ can be understood as a systematic fractional instanton expan- 
sion, where the power of Si determines the number of fractional instantons in the factor 
U(Ni) we are taking into account. In this note we will be considering just the one instanton 
sector, which is determined by the universal term coming from the free (conifold) theory. 
We leave the study of multi-instanton effects due to the J-q expansion for future work. 

2.2 D-brane instantons 

Let us now proceed to review the aspects of D-brane instanton calculus most relevant to 
our discussion. As reviewed in the previous section, the Calabi-Yau has cycles over which 
we can wrap space-filling D-branes in order to engineer 4d supersymmetric gauge theories. 



5 When computing instanton effects, sometimes there will be no spacetime filling D-brane wrapping the 
cycle, and thus the definition of r as a gauge coupling constant does not make sense. The definition in 
terms of the (complexified) volume of the resolved cycle still makes sense, and it is the one we will always 
be implicitly using. 
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Depending on the particular vacuum we choose there might also be cycles not wrapped 
by any spacetime filling D-brane. Over these cycles we might wrap D-brane instantons, 
which in this context are branes wrapping these internal cycles and localized at a point in 
spacetime. They give nonperturbative (in g s ) contributions to the 4d effective action living 
on the spacetime wrapping branes. 

The basic condition for a D-brane instanton to contribute to the superpotential in 
this setting is that it has exactly two unlifted fermionic zero modes, which will play the 
role of the 9 variables in the integration of the superpotential over superspace. The kind 
of backgrounds where are considering preserve, in the absence of branes and orientifolds, 
8 supercharges. A D-brane instanton is a 1/2 BPS object, and thus has 4 fermionic zero 
modes 9 a , 9a corresponding to the 4 spontaneously broken supersymmetries. If the D-brane 
instanton is going to contribute to the superpotential, we should include some mechanism 
that lifts the two 9 a Goldstinos. 

The simplest configurations in which this happens, and which have no simple gauge 
theory interpretation, are instantons of the U(l) and Sp(0) types. We describe them in 
the following. 

An instanton of Sp(0) type occupies a cycle in which there are no wrapped space-filling 
branes, and which is mapped to itself by the orientifold action. Furthermore, the sign of 
the orientifold projection is such that the gauge group on the instanton is 0(1). Such a 
projection does indeed remove the 9 a Goldstinos, as we want 6 . 

In order to make the connection with the matrix model more transparent, here we have 
chosen to refer to this kind of instanton as u Sp(0)", instead of using the more common 
denomination "0(1)". The justification for this denomination from the instanton point of 
view is the following. Consider the orientifold action on N space-filling D-branes wrap- 
ping the same cycle as the instanton. Since the orientifold projection acts oppositely on 
space-filling and instanton branes, these branes would have gauge group Sp(N). The D- 
brane instanton in this context can be interpreted as the gauge instanton for Sp(N), with 
the couplings for the string modes between the different branes implementing the ADHM 
construction [28, 29, 30, [32]| . The notation Sp(0) is thus intended to suggest that the 



D-brane instanton can be thought of as a gauge instanton of a (admittedly rather empty) 
Sp(0) gauge group. 

Let us now describe instantons of the U(l) type. These instantons do not exist in 
pure field theory either, but they do exist in string theory when there is a single D-brane 



wrapping the same cycle as the instanton [14, 11, 33 1. 



In order to explain why this happens, let us consider first N spacetime filling branes 
wrapping the cycle, giving a U(N) gauge theory. Intuitively, since it is on top of the D- 
branes, the D-brane instanton "feels" just J\f = 1, and thus it has the two desired Goldstinos 
9 a only. More precisely, the couplings between the zero modes of the instanton and the 
fields living in the branes implement the ADHM constraints, and the ADHM couplings lift 
the 9 a modes. 



6 See for example |?7] for a nice discussion of the details of the possible orientifold actions and their effect 
on the instanton zero modes 
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One interesting observation one could make is that the way in which the couplings 
between the zero modes implement the ADHM constraints and lift the two extra goldstinos 
does not require that N is bigger than 1 for U(N), it works for U(l) gauge theories too 
as studied in detail in [33]. U(l) theories do not have instanton effects in gauge theory, 



and the fact that D-brane instantons contribute in this case is very much related to the 
fact that we are completing this £7(1) theory in the UV by embedding into string theory. 
To restate one of the main points of || and also of this work: such a completion is not 
as innocuous as it looks, since it also changes the low energy superpotential due to these 
exotic U(l) instantons. 

Before going on to check the agreement between the matrix model and the D-brane 
instanton calculus in these cases, let us make a couple of important remarks. Although we 
restrict here to the case in which the instanton has exactly two zero modes, the general 
case in which the instanton has more zero modes is also interesting. What happens here is 
that the instanton does not contribute to the superpotential, but can generate couplings 
involving more than two fermions, called higher F-terms [|34| , |3q| . One important example 
is the gauge instanton for SQCD with Nf = N c . These instantons do not contribute to 
the superpotentials, but induce operators with insertions of more than two fermions (or 
equivalently, derivatives). They are of the form: 

J d i xd 2 ew- i - h ___ In j n ($, $) (pW 1 ■ W 1 ) • • • (pW n ■ W") (2.16) 

where <3> represents the moduli of the system, and w is a section of f?^ ® ^?vf, with A4 
the moduli space 7 . We include some preliminary remarks on the relation between matrix 
models and higher F-terms in the conclusions. 

Another interesting possibility for lifting extra zero modes which we have not discussed 



here is by using fluxes, as studied in pa, 37, 38, 3£, M. Modifying the matrix model side 



of the story in order to incorporate the effect of these fluxes is a worthwhile and interesting 
problem which we will not attempt to solve here. 

3. The glueball superpotential for low rank and D-brane instantons 

In this section we will study the two basic cases in which the matrix model result does 



not match the gauge theory result, namely, as reviewed in Section 2.1, configurations with 
gauge group U(l) and Sp(0). We will see that there is a nice agreement between the 
result of the D-brane instanton calculation and the result obtained from the matrix model 
prescription. This gives evidence for the advertised result that matrix models can be used 
to compute D-brane instanton effects. 

The geometry we will be focusing on will be the resolved conifold, so let us take a 
quadratic superpotential for the gauge theory: 

d 2 6TrW{<$>) = f d 2 8Tr^- (3.1) 



7 In other words, it is a form antisymmetric in the i indices, antisymmetric in the j indices, and these 
indices live in the cotangent bundle over the moduli space. 
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where the factor of 1/2 is for later convenience. Plugging this W(x) into equation ( |2.7| ) 
we obtain the conifold equation x 2 + y 2 + z 2 + w 2 = 0. As explained above, we will be 
considering the resolved conifold, with a resolution parameter determined by the theory 
we are engineering. 

The matrix model computation 

Let us start by the matrix model side of the story, which is the simplest. We have to 
compute the planar limit of the matrix model free energy in this background. The conifold 
is described by the free holomorphic one-matrix model, so we just obtain the Veneziano- 
Yankielowicz result. Note that we must set iV = 1 in eq. (|2.5|) in order to capture Sp(0) 



and U(l) D-brane instanton effects, as discussed in Section 2.1. We get: 



W eff (S, t)=S log (^pj + S - 27TirS (3.2) 

with Ao the cutoff scale and r the bare coupling of the configuration, defined in eq. (|2.9| ): 

2vrir = — / J + iB (3.3) 
9s Js2 

We can now proceed to integrate S out. Its equation of motion gives: 

S = e~ 2mT Al (3.4) 
which upon substitution into ( |3.2| ) gives: 

W = e- 2niT A 3 = A 3 (3.5) 
where we have introduced the cutoff independent scale A. 
The D-brane instanton computation 

Let us now describe in some detail how the stringy instanton result reproduces this con- 
tribution. In order to avoid the complications of having to impose the (super)ADHM 
constraints, we will focus on the Sp(0) case, so the 9 modes are simply projected out by 
the orientifold. A good place to see how the saturation of the 9 modes goes for U(l) 
instantons is for example |33|| . 

We will start by considering the case with W = 0. In this case we will have 8 back- 
ground supercharges close to the orientifold, and 16 in the bulk. The Sp(0) instanton is 
a 1/2 BPS object wrapping a 2-cycle mapped to itself under the orientifold action, so it 
will have 4 Goldstinos in its worldvolume. Let us call the surviving Goldstinos 9^, where 
A is an index taking values in the fundamental representation of the SU(2) R-symmetry 
group preserved by the orientifold. Such an instanton has four neutral zero modes, and 
thus cannot contribute to the superpotential. 

We now switch on the quadratic superpotential W(z) = \mz 2 . This reduces the 
supersymmetry of the background by half, so now we have 8 supercharges on the bulk and 
4 close to the orientifold. In such a configuration it is reasonable to guess that two of the 
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four zero Goldstinos of the instanton get lifted, and the instanton can now contribute to 
the superpotential. This is indeed the case, let us now explain how this deformation acts 
from the point of view of the effective action on the instanton zero modes. 

Let us start by studying the case were we have background spacetime filling branes, so 



we are dealing with ordinary gauge theory instantons. As reviewed in Section 2.1 , in the 
gauge theory the deformation of the background appears as a breaking of the supersym- 
metry from Af = 2 to Af = 1 by the addition of the superpotential for the adjoint chiral 
multiplet: 

W($) = imTr$ (3.6) 

One consequence of such a coupling is that some of the zero modes associated to the field 
<3? get lifted, let us call such modes \ a and Aq,. We have chosen this notation in order to 
remind the reader that the adjoint superfield comes from the vector multiplet of Af = 2, so 
the highest component fermion of the multiplet can be thought of as a gaugino from the 
point of view of the broken Af = 1 algebra. As reviewed in Section VI. 4 of [41|, this can be 



nicely encoded in the effective action for the instanton zero modes as the addition of the 
following term: 

2 

fiSi ns t = A"Aq, (3-7) 

9YM 

Such a term can be used to saturate the d 2 X term in the instanton measure, allowing the 
instanton to contribute to the superpotential. 

We now claim that a term similar to ( |3.7D also exists for our Sp(0) instanton. We 
cannot compute the existence of such a coupling via CFT methods in string theory, but 
given the fact that string theory gives rise to the whole ADHM construction in such a 



natural way [42, 43], with a natural identification of the instanton- instanton strings as 
gauge theory zero modes, assuming that an analog of ( fTjj ) holds seems to be reasonable. 

From the point of view of the instanton brane X a is one of the Goldstinos 9£ , let us 
take it to be Q 2 a . The coupling induced by the deformation would then be: 



-— fl 2 ^ (3.8) 

9YM 



Just as in the gauge theory case, we can bring down this coupling in order to saturate the 
fermionic integration over the Q\ modes. This integration will give us a factor of m in front 
of the resulting contribution to the superpotential (the factors of it and qym get taken care 
of once we integrate over the rest of the zero modes properly |ftlfl). 

The rest of the calculation, once we have dealt with the massive adjoint, is by now 
well understood (see for example f33[l for a careful derivation for the more involved U(\f) 
instanton). The end result is simply given by: 

W ef f(T) = mMl exp (-2vrir) (3.9) 



The M s factors come from the measure of integration of the zero modes [£g, [3^, 33], and 



the exponential suppression comes from the volume of the cycle wrapped by the instanton. 
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Figure 1: Quiver for the orientifolded orbifold of the conifold we discuss in the main text. We have 
omitted the part of the quiver to the right of U(N2), which will not be relevant for our discussion. 

In order to match with the matrix model result we need to relate the low energy 
physical A appearing in ( |3.5D with the microscopic scales appearing in (|3 . 9|) . This goes 
as follows. The physical scale for the theory before integrating out the adjoint is given in 
terms of the string scale quantities as: 

A 2 adj = M 2 e~ 2 ™ (3.10) 

The m prefactor gives the matching relation between the theory before and after integrating 
out the adjoint: 

A 3 = A 2 adj m (3.11) 
The resulting superpotential can thus be written as: 

W = A 3 (3.12) 

matching (3.5) beautifully. 



4. A quiver example: Z n orbifold of the conifold 

Let us proceed to the more interesting case in which the D-brane instanton intersects some 
space-time filling branes. In this discussed in Jj], ||, g|, there are fermionic zero 

modes arising from strings stretching between the instanton and the brane. Integration over 
them typically induces interesting couplings in the world volume theory of the brane, which 
might be potentially relevant for model building. They generically give rise to couplings 
with are forbidden in perturbation theory. 

In this section we will focus on a model which is easy to analyze using matrix model 
techniques. It is the theory corresponding to branes at the singularity of an abelian orbifold 
of the conifold. These theories can be easily obtained from A n orbifold quivers by giving 
large masses to the adjoints. The analysis we do here is easily generalizable to Af = 2 quiver 
theories with superpotentials for the adjoints, as in the original paper of Dijkgraaf and Vafa 
@. We will choose an assignment of ranks of the quiver such that the orientifolded node 
is empty, giving rise to one of the Sp(0) factors we have been describing. 

The relevant quiver diagram is shown in Figure |l|, where we are omitting the part of 
the quiver that will not be relevant for us. The nodes of the diagram denote gauge factors, 
and the arrows denote bifundamental chiral multiplets. The theory also has a quartic 
superpotential, given for generic ranks in the nodes by: 

W = TiX 01 X 12 X 21 X 10 (4.1) 
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In our case, since we are setting the rank of the Sp node to zero, Xoi and X\q vanish, 
and W = 0. Nevertheless, as we saw in Section \2,.1[ the ranks that go into the matrix 
model computation are unrelated to those in the physical string. This means that we can 
compute the effective superpotential keeping Xqi, X\q and W in the game, and treating 
them as if we had Sp{K) with K > 0. The only place where the actual rank of Sp{K) 
enters in the matrix model computation is in determining the prefactor in the Veneziano 
Yankielowicz term, as in eq. (|2.5|). 



Let us proceed to compute the low energy superpotential from the matrix model point 
of view. We want to obtain the effect induced by one instanton in the Sp(0) node. One 
way we could go about computing this effect would be by applying the prescription in Q 
for A n quivers and taking the limit where the adjoint masses go to infinity, keeping the 
dynamical scales of the resulting conifold theory finite. See for example |44], ^5], |46|, ^7]] for 
calculations along these lines. Since we are not interested in the full answer, but just the 
leading term due to Sp{K), we can take instead the following shortcut. 

The leading term we are interested in corresponds to letting Sp{K) confine while taking 
the rest of the quiver theory to be an spectator. In this way XJ(N\) x U(N 2 ) can be taken 
to be a flavor group, and X\2, X21 are background vevs. From this point of view, we can 



take (4.1) to mean: 

W = TiX m X l2 {X 2l X l2 ) (4.2) 

In this equation {X 2 iX\ 2 ) is the mass matrix for Xq\, -X12, which we can now integrate 
out. It is important to take into account the matching relation: 

Aff +1) = Ajg* 1 '-" det(X 21 X 12 ) (4.3) 

where Khigh and are the dynamical scales of the Sp{K) factor before and after inte- 
grating out the massive matter. 

The remaining theory is Sp{K) with no flavors. It confines, and the low energy super- 
potential is given by: 

W = (K + 1)AL = (K+ l)A 3 ^f /{A ' +1) det(X 21 X 12 )^ (4.4) 

We can now use this result computed in the gauge theory with K > 0, extract the 
matrix model result for dZo/dS, and apply it to the Sp(0) case. We obtain the result: 

W eff = A 3 "^ det{X 21 X 12 ) (4.5) 

Let us now compare this result with the D-brane instanton computation. The tech- 
nology for dealing with this problem is by now standard |7], |8|, |9[. An instanton sitting 
on the Sp(0) node has two neutral zero modes 8 a , and two zero modes a, j3 (these are 
anticommuting scalars living in the fundamental and antifundamental of U{N\)) arising 
from strings going from the instanton to the U{N\) brane. Essentially the same physics 



1 48] that gives rise to the superpotential ( [4.1[) gives rise to the following coupling in the 
instanton action: 

5S inst = aX 12 X 2l p (4.6) 
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The contribution of the instanton to the low energy effective superpotential is then 
given by: 

W eff = mM 2 - N ^e- 2mT J dad/3 exp (5S inst ) (4.7) 

The overall prefactor works exactly as in the conifold case studied in Section |3|, and we can 
identify it with A 3 ~ Nl above. The second term can be used to saturate the integrations 
over a and (3, and gives rise to: 



/ 



dad/3 exp (8S inst ) = detX 12 X 2 i (4. 



reproducing the result (4J5) obtained from the matrix model. 



We see that the essential ingredients that went into the previous calculation are those 
that appeared already in Section ||, where we studied the case of the isolated instanton. 
Namely, once we considered the glueball superfield for the node containing the instanton 
dynamical, the matrix model correctly gave us the same answer as the (in general more 
involved) instanton calculation. 

5. Conclusions and further directions 

In this work we have started the study of D-brane instanton effects from the point of view of 
matrix models. We found that at least the simplest examples indicate that, in the same way 
that ordinary gauge non-perturbative effects are captured by a perturbative calculation in 
the matrix model, more stringy effects are also computed by the matrix model. We checked 
this result in some detail for isolated conifolds, and for quiver gauge theories admitting a 
matrix model description. 

There are many ways in which one could extend the results in this note, let us list here 
a few of them. 

We have restricted ourselves to matching the matrix model and the D-brane instanton 
results just up to leading order in the instanton expansion, namely the one (fractional) 
instanton sector. This simplified the calculations enormously, reducing the problem to 
writing down the Veneziano-Yankielowicz superpotential for the instanton node and ex- 
tremizing the complete superpotential. This is already a good sign, and indicates that the 
matrix model formalism might be significantly simpler to deal with (when it is applicable) 
than the D-brane instanton calculus, which is sometimes involved even in the one-instanton 
sector. 

In this same line of argument, one technical difficulty to match the higher orders in 
the matrix model perturbation expansion with non-perturbative effects in string theory is 
that we will have to deal with multi- (fractional) instanton effects in the string side, and 
these are somewhat cumbersome to deal with using standard (D-brane) instanton calculus. 



See however [41] for a review of some powerful techniques for dealing with multi-instantons 



in the case of gauge theory with extended susy, |4{J for progress in the M = 1 case, and 
[36, 11, |50|, 51, [37], [5^, 53] for some recent progress in more stringy systems. It would 
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be interesting to check the proposed correspondence for higher orders in the instanton 
expansion using some of these techniques. 

The geometries on which one can apply directly the ideas of Dijkgraaf and Vafa are 
rather restricted, and not very suitable for model building 8 . Since the time when the DV 
results appeared there have been very interesting developments in the study of matrix 
models and topological string theory. Particularly relevant for our story is ^4|, where a 
matrix-model inspired formalism is given for computing perturbative amplitudes of the 
B- model on a very rich set of geometric backgrounds (in particular, this includes mirrors 
of toric Calabi-Yau spaces). It would be extremely interesting to see if and how D-brane 
instanton effects appear in this setup. 

We would also like to make some preliminary comments on the study of higher F- 
terms in the context of matrix models (we have succinctly introduced these higher F-terms 



towards the end of Section 2J2). Recently it has become clear that these higher F-terms 
form an integral part of the global picture for D-brane instanton effects across moduli space 
1 11, 40]. They have a fascinating interplay with ordinary superpotential terms in order 



to give a consistent description of non-perturbative effects across the whole moduli space of 
the compactification. It would be rather interesting to see if and how these higher F-terms 
are encoded in the matrix model description. 

One seemingly reasonable possibility is that these terms are computed by the non- 
planar diagrams in the matrix model. A naive and straightforward reason is that we need to 
compute terms in the effective superpotential with derivatives, or more than two fermions, 
and it is well known that in the topological string higher genus diagrams are associated with 
operators of this kind in the physical superstring. This is obviously too sketchy, so it would 
be good to have a more direct connection between non-planar diagrams and Beasley-Witten 
higher F-terms. We know that non-planar diagrams encode the effect of the C-deformation 



of the gluino anticommutation algebra [55|, so if non-planar diagrams are also related to 
higher F-terms then we would have a nice relation between Beasley-Witten higher F-terms 
and the C-deformation. 

As a last remark, all of our discussion has been in the context of type IIB string 
theory. By mirror symmetry, one could expect that Chern-Simons theory in the appropriate 
background might also compute the effect of exotic E2 instantons. It would be rather 
interesting to see how this works. 
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